Difference in nature of correlation between 
NASDAQ and BSE indices 



P. Manimaran 

School of Physics, University of Hyderabad, Hyderabad 500 046, India 

Prasanta K. Panigrahi 1 ' 2 , and Jitendra C. Parikh 2 

1 Indian Institute of Science Education and Research (Kolkata), Salt Lake, Kolkata 

700 106, India 

2 Physical Research Laboratory, Navrangpura, Ahmedabad 380 009, India 
Abstract 

We apply a recently developed wavelet based approach to characterize the corre- 
lation and scaling properties of non-stationary financial time series. This approach 
is local in nature and it makes use of wavelets from the Daubechies family for de- 
trending purpose. The built-in variable windows in wavelet transform makes this 
procedure well suited for the non-stationary data. We analyze daily price of NAS- 
DAQ composite index for a period of 20 years, and BSE sensex index, over a period 
of 15 years. It is found that the long-range correlation, as well as fractal behavior 
for both the stock index values differ from each other significantly. Strong non- 
statistical long-range correlation is observed in BSE index, whose removal revealed 
a Gaussian random noise character for the corresponding fluctuation. The NASDAQ 
index, on the other hand, showed a multifractal behavior with long-range statistical 
correlation. 
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1 Introduction 



Non-stationary time series have been investigated through several approaches. 
In particular, the characterization of fluctuations and their scaling behavior 
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have been the focus of many studies, since they reveal the nature of the dynam- 
ics. In this context, financial time series have attracted considerable attention 
(Q]). The large length of the available data of various stock market indices 
make them ideal candidates for analysis. Furthermore, the complex dynamics 
of the variations in stock prices yield fluctuations which can show correlations, 
as well as scaling behavior. The goal of the present paper is to study the self- 
similar and correlation properties of NASDAQ and BSE indices belonging to 
two different economical environments. These stock indices belonging to de- 
veloped and developing countries may show characteristic differences, possibly 
arising due to differences in their underlying dynamics. We concentrate on the 
nature of correlations and fractal behavior for which wavelet transform (0; |i|) 
is used as a tool for extracting the fluctuations at different scales. 



A number of methods have been devised to find scaling behavior in time se- 
ries. The well-known structure function method (0) and the recently developed 
wavelet transform modulus maxima (WTMM) method (js)), relying on contin- 
uous wavelet transforms, are widely used for the analysis of stationary data. 
The fact that most of the time series arising in real systems are non-stationary 
in nature introduces complications in estimating the scaling behavior, while 
using the above two approaches, which are global in nature. Hence, in re- 
cent times, local approaches, like detrended fluctuation analysis (DFA) and 
its generalization MF-DFA (jj; 0; @; 0; QiJ) have been developed to handle non- 
stationary data. In this case, one uses windows of various sizes to separate 
fluctuations from the trend, which can also be shuffled to remove any corre- 
lation in the data set. While isolating the average or the trend of the data 
points in a given window, one takes recourse to linear or quadratic fit in the 
DFA a ppro ach. We have introduced a new method based on discrete wavelets 
( 111 ; Il2j ; Il3l ) to characterize the scaling behavior of non-stationary time series. 



The present procedure is similar to those in MF-DFA (lid ), except that in 
order to detrend, we use wavelets and MF-DFA uses local polynomial fits. 
Recently, this method has been used by Brodu, to analyze in real time, fractal 
behavior of dynamic time series (11 4 ). The relative merits of MF-DFA and a 
variety of other approaches to characterize fluctuations have been carried out 
in Ref . (flil) . It is worth emphasizing that fluctuation analysis and character- 
ization have been earlier attempted using Haar wavelets, in the context of 
bio-medical applications, without the study of scaling behavior (fl6l ; Il7l ). 



Wavelets from Daubechies family are used for extracting trend from the given 
data set. Fluctuations are captured by high-pass coefficients and the trend cap- 
tured by the low-pass coefficients of wavelet transform. The discrete wavelet 
transform provide a handy tool for isolating the trend in a non-stationary data 
set, because of its built-in ability to analyze data in variable window sizes. In 
this note, we analyze returns of stock index values through our new wavelet 
based method. Multi-fractal properties are also investigated using multi-fractal 
spectrum. We analyze daily price of NASDAQ composite index for a period 



2 




500 1000 1500 2000 2500 3000 3500 4000 4500 5000 




500 1000 1500 2000 2500 3000 3500 4000 4500 5000 



r 1 1 n i 

















500 1000 1500 2000 2500 3000 3500 4000 4500 5000 
time (days) 



Fig. 1. [a] NASDAQ daily (close) composite index for a period of 20 years, starting 
from ll-Oct-1984 to 24-Nov-2004, [b] daily returns show more clusters of small and 
large fluctuations and [c] the returns after shuffling show disappearance of clustering 
behavior. 

of 20 years, starting from ll-Oct-1984 to 24-Nov-2004, and BSE sensex index, 
over a period of 15 years, starting from 2-Jan-1991 to 12-May-2005. 



2 Data analysis 



It is found that the nature of correlation is quite different between these two 
financial time series and significant non-statistical correlation exists in both of 
them. Removal of the same reveals that BSE index is primarily mono-fractal 
with the fluctuations showing a Gaussian random noise character. On the 
other hand, the NASDAQ index shows a weak multifractal behavior with long 
range statistical correlation. 

From the financial (NASDAQ composite index and BSE sensex) time series 
x(t), we first compute the scaled returns defined as, 

G(t) = [ln(x(t + l))-ln(x(t))]/a, t = 1, 2...(N - 1); (1) 

here a is the standard deviation of x(t). From the returns, the signal profile 
is estimated as the cumulative, 

Y(i) = ip(t)], i = l,....,N-l. (2) 

t=i 

Next, we carry out wavelet transform on the profile Y(i) to separate the fluc- 
tuations from the trend by considering precise values of window sizes W cor- 
responding to different levels of wavelet decomposition. We obtain the trend 
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Fig. 2. [a] BSE sensex daily (close) index value for a period of 15 years, starting from 
2-Jan-1991 to 12-May-2005, [b] daily returns show much less (compared to Fig. 1 
(b)) clustering of fluctuations and [c] the returns after shuffling show no significant 
difference in appearance relative to Fig. 2 (b). 



by discarding the high-pass coefficients and reconstructing the trend using 
inverse wavelet transform. The fluctuations are then extracted at each level 
by subtracting the obtained time series from the original data. Though the 
Daubechies wavelets extract the fluctuations nicely, its asymmetric nature and 
wrap around problem affects the precision of the values. This is corrected by 
applying wavelet transform to the reverse profile, to extract a new set of fluc- 
tuations. These fluctuations are then reversed and averaged over the earlier 
obtained fluctuations. These are the fluctuations (at a particular level), which 
we consider for analysis. In Figs. 1 and 2, we give the time series for the two 
index data sets, and the corresponding returns. We also show shuffled returns 
for the two series to examine the correlation as well as "bursty" (clustering) 
behaviour. 

The extracted fluctuations are subdivided into non-overlapping segments M s = 
int(N/s) where s = 2 ( - L ~ 1 ^W is the wavelet window size at a particular level 
L for the chosen wavelet. Here W is the number of filter coefficients of the 
discrete wavelet transform basis under consideration. For example, with Db-4 
wavelets, s = 4 at level 1 and s = 8 at level 2 and so on. It is obvious that 
some data points would have to be discarded, in case N/s is not an integer. 
This causes statistical errors in calculating the local variance. In such cases, 
we have to repeat the above procedure starting from the end and going to 
the beginning to calculate the local variance. The detrending and extracted 
fluctuations have been depicted in Figs 3 and 4. 

The q th order fluctuation function F q (s) is obtained by squaring and averaging 
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Fig. 3. [a] Detrending the integrated returns of NASDAQ composite index at the 
scale level-3, through Db-8, wavelet, and [b] the extracted fluctuations at the scale 
level-3 (window size 32). 
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Fig. 4. [a] Detrending the integrated returns of NASDAQ composite index at the 
scale level-4, through Db-8, wavelet, and [b] the extracted fluctuations at the scale 
level-4 (window size 64). 

fluctuations over all segments: 

i 2M S 

F q {s) = {—Y.[F 2 (b,s)}^Yl<i. (3) 

ZIV1 s b=l 



Here 'g' is the order of moments that 
is repeated for variable window sizes 
The scaling behaviour is obtained by 

F q (s) ~ s h ®, 



takes real values. The above procedure 
for different values of q (except q — 0). 
analyzing the fluctuation function, 

(4) 



in a logarithmic scale for each value of q. If the order q = 0, direct evaluation 



5 



0.8 

0.7 
_ 0.6 
^ 0.5 

0.4 

0.3 

-10 -8 -6 -4 -2 2 4 6 8 10 

q 

6 
4 
2 

a 
*" -2 
-4 

-6[ 

-8 C 

-10 -8 -6 -4 -2 2 4 6 8 10 

q 

Fig. 5. [a] (NASDAQ composite index) Scaling exponents h(q) values for various 
q values and [b] r(q) representation of h(q) values for various q values, where 
r{q) = qh(q). 

through Eq. (3) leads to divergence of the scaling exponent. In that case, 
logarithmic averaging has to be employed to find the fluctuation function: 

F q (s) = exp{-±- £ ln[F\b, s)]" 2 } 1 ". (5) 

ZJVi s 6=1 



3 Results and Discussion 

As is well-known, if the time series is monofractal, the h(q) values are indepen- 
dent of q. For multifractal time series the h(q) values depend on q. The cor- 
relation behaviour is characterized from the Hurst exponent (H = h(q = 2)), 
which varies from < H < 1. For long range correlation, H > 0.5, H = 0.5 
for uncorrelated and H < 0.5 for long range anti-correlated time series. 

The scaling exponent is calculated for various values of q for both the stock 
index values. Figs. 5 and 6, show the way h(q) and r(q) vary with q for the 
returns and shuffled returns of the two time series. The non-linear behaviour 
of h(q) for different q values, is the measure of multifractality. 

The scaling behaviour of the observed data sets can also be studied by evaluat- 
ing f(a) spectrum. f(a) values are obtained from Legendre transform of r(q): 
f(a) = qa — r(q), where a = For monofractal time series, a = const., 

whereas for multifractal time series there occurs a distribution of a values. 
The f{ct) spectra for the two time series are shown in Figs. 7 and 8. For the 
unshuffled returns, one observes a broader spectrum, whereas for the shuffled 
returns, where the correlation is lost, the same is narrower. 
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Fig. 6. [a] (BSE sensex) Scaling exponents h{q) values for various q values and [b] 
r{q) representation of h(q) values for various q values, where r(q) = qh(q). 




Fig. 7. From the integrated returns of NASDAQ composite index values, the calcu- 
lated multifractal spectrum is broader than the spectrum of shuffled returns. 




Fig. 8. From the integrated returns of BSE sensex values, the calculated multifractal 
spectrum is much broader than the spectrum of shuffled returns, when compared 
to the NASDAQ composite index values. 
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Fig. 9. The semi- log plot of distribution of logarithmic returns of NASDAQ com- 
posite index values compared with Gaussian distribution. 

The semi-log plot of distribution of logarithmic returns of NASDAQ composite 
index values is shown in Fig. 9. It exhibits fat tails and non-Gaussian features. 
In case of BSE sensex, the semi-log plot of distribution of logarithmic returns 
shows fat tails, which are less prominent, as shown in Fig. 10. The distribution 
is quite similar to Gaussian white noise, revealing distinct differences between 
NASDAQ composite index and BSE sensex values. Although correlation is 
present in the two time series, they reveal distinct probability distributions 
once the correlation is removed. Furthermore, we observer that in case of BSE 
sensex values, the calculated multifractal spectrum is much broader than the 
spectrum of shuffled returns, when compared to the NASDAQ composite index 
values. 



4 Conclusion 



In conclusion, the wavelet based method presented here is found to be quite 
efficient in extracting fluctuations from trend. It reveals the distinct differences 
in the long-range correlation, as well as fractal behavior of the two stock index 
values. Strong non-statistical correlation is observed in BSE index, whereas 
the NASDAQ index showed a multifractal behavior with long-range statis- 
tical correlations. In case of BSE index the removal of correlation revealed 
the Gaussian random noise character of the fluctuations. It is interesting to 
note that the effect of country specific parameters like corruption on economic 
development and investment has been recently quantified through scaling anal- 
ysis (1181) . In a similar manner, the above observed differences between the two 
stock indices belonging to two different economic environment is probably due 
to local dynamics. 
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Fig. 10. The semi-log plot of distribution of logarithmic returns BSE sensex values 
compared with Gaussian distribution. 
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